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The electromagnetic form factors of the A baryon are evaluated within the framework of a co- 
variant spectator quark model, where S and D-states are included in the A wave function. We 
predict all the four A multipole form factors: the electric charge Geo, the magnetic dipole Gmi, the 
electric quadrupole Ge2 and the magnetic octupole Gm3- We compare our predictions with other 
theoretical calculations. Our results are compatible with the available experimental data and recent 
lattice QCD data. 



I. INTRODUCTION 

In the history of hadronic physics the nucleon was one 
of the first particles to have its internal structure exper- 
imentally uncovered. The nonpointlike character of the 
nucleon was suggested by the measurements of the pro- 
ton anomalous magnetic moment. Around 1950, mea- 
surements disclosed an exponential falloff for the charge 
distribution, and three decades later the SLAC measure- 
ments showed that the corresponding charge form factor 
in momentum space coincided, almost perfectly, with the 
magnetic dipole form factor. This image of the proton 
proved to be a good approximation for low Q 2 , the square 
of the four-momentum transfer, and was only superseded 
by the results from the Jlab polarization measurement in 
1999 [|. 

The study of the next baryonic state (the A) has been 
more challenging. Although the pure spin 3/2 and isospin 
3/2 structure of the A was clearly isolated from the back- 
ground of the pion nucleon scattering experimental cross 
sections, its theoretical description is not yet totally clear. 
Quantum chromodynamics (QCD) emphasizes the three- 
quark state component, and historically the existence of 
the A ++ led to the introduction of color as a quantum 
number. Still, the issue on how much of the A struc- 
ture comes from a three-quark state, and how much is a 
molecularlike state of a nucleon and a pion, remains an 
open question. 

Since the A is a spin 3/2 particle, its electromag- 
netic structure can be characterized by four form factors, 
namely the electric charge Geo, magnetic dipole Gmi, 
the electric quadrupole Ge2 and electric octupole Gmi 
form factors. The first two describe the charge and mag- 
netic dipole distributions, while the last two measure the 
deviation of the first ones from the symmetric form (2j. 
At zero transferred momentum squared, Q 2 — 0, the form 
factors define the A multipole moments: the charge eA, 
the magnetic moment /iA, the electric quadrupole mo- 



ment Qa and the magnetic octupole moment Oa • If the 
A would simply be a completely symmetric state, with 
no quark D-states relatively to the quark pair (diquark), 
both Qa and Oa would vanish 

The unstable character of the A (with a mean lifetime 
of 5.6xl0 _24 s) makes difficult to probe the electromag- 
netic properties of this particle/resonance. Until recently 
the available information on the A structure came from 
the study of the 7-/V — >• A reaction only. The dominance 
of the magnetic dipole moment in that reaction asserts 
the spin flip of a quark in the nucleon as the main pro- 
cess, but the magnitudes of the quadrupole contributions 
(G* E and Gq) indicate a small deformation of the A Q. 

The experimental information on the A is restricted 
to the A ++ and A + magnetic moments (5l-fl2j. Even so, 
those results are affected by considerable error bars, due 
to the unstable nature of the A and theoretical model 
uncertainties. This is particularl y tr ue for the A + mag- 
netic moment Q. See Refs. @, H3 for a review. New 
experiments are in progress at MAMI [l2| to extract 
/iA+ more accurately with the help of reaction models 
based on hadronic degrees of freedom [ll|, 0, E3 ■ The 
scarce experimental information on /^a is compensated 
by extensive theoretical studies within various frame- 
works [ll|, 14 50]. For a review see Ref. • As for other 
A observables, there is no other information, with the 
exception of the electric quadrupole moment Qa , which 
was estimated usin g a connection to the jN — > A tran- 
sition quadrupoles [5l|. For Q 2 ^ there is no experi- 
mental information related to the A form factors. Even 
a very basic property like the A electric charge radius 
is at present rather unknown, contrarily to what hap- 
pens for the nucleon, where that quantity is accurately 
measured [52]. Our knowledge on the A charge radius 
relies so far only on diverse theoretical model calcula- 
tions [H, iH, M, HI El [H-i2. There are also the- 
oretical calculations of the A quadrupole moment Qa 
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Recently, the nature of the A was explored in another 
direction, and under another light. Lattice QCD, the dis- 
crete version of the fundamental nonperturbative theory 
of the hadrons, was used to estimate the A elastic form 
factors [731 - 176] ] , following the pioneer works of Bernard et 
al. [z3 and Leinweber et al. The works [zl[7l|7| 
are based on the quenched approximation. The works 
[74l [75| are unquenched calculations, which include also 
sea quark effects. The results of the MIT-Nicosia group 
[73M75I were obtained for a wide r ang e of Q 2 , while the 
calculations of the Adelaide group [76[ are performed for 
only one small value of Q 2 (the limit case Q 2 = is 
impossible to reach for technical reasons) and extrapo- 
lated for Q 2 — to obtain the A multipoles and the 
charge and magnetic radii. In particular for Q 2 = 
the background field method was used to calculate with 
a good accuracy the A magnetic moment [ij [zH l80| - 
These methods were applied to pion mass values within 
the 0.3—1 GeV region and consequently some extrapo- 
lations are required 0, [8l|. In the absence of direct 
experimental information lattice QCD provides the best 
reference for theoretical calculations for finite Q 2 . 

Motivated by the recent publication of lattice QCD 
data for the A form factors p73MT6l |80T | , several models 
were also extended to determine those physical quanti- 
ties. Calculations of the A form factors were performed 
within the covariant spectator theory 0, [HJ and also 
within a chiral Quark-Soliton model (%QSM) [62]. Also, 
the magnetic octupole moment 0\ was estimated by 
Buchmann [2j, using a deformed pion cloud model, and 
by Ceng (49J, using xPT. Both Qa and Oa were esti- 
mated using QCD sum rules [69| and a covariant specta- 
tor quark model [82j |. 

In this work we use a constituent quark model based 
on the covariant spectator formalism [4, 83 85] to predict 
the Q 2 dependence of all four A multipole form factors. 
Our results will be compared with the available exper- 
imental data, and the results from other models in the 
literature, as well as with recent lattice data. The A 
electromagnetic form factors give us an insight into the 
A internal structure. In particular the determination of 
the A electric quadrupole and magnetic octupole mo- 
ments gives information on how the A is deformed, and 
on its size and shape. And this information is valuable, as 
we will conclude, to adequately constrain quark models. 
These in turn can be used to extrapolate and interpret 
lattice results. 

The nucleon, as a spin 1/2 particle, has no electrical 
quadrupole moment, and the A, a spin 3/2 state, emerges 
as the first candidate for a deformed baryon. A deforma- 
tion (compression or expansion along the direction of the 
spin projection) would imply Qa ^ 0. The interpreta- 
tion of Qa ^Oasa signature of distortion matches well 
with our intuitive notion of deformation, based on the 
nonrelativistic limit of the Breit frame form factors. In 
that limit the Fourier transformation of the form factors 
gives the spatial distribution of charge or magnetic mo- 
ment. We will therefore identify this notion as classical 



deformation. Recently, an alternative interpretation of 
distortion, based on transverse densities defined in the 
infinite moment frame, was introduced by other authors 
[75l [86- 88] . Our calculations in this paper provide the 
basis needed for a forthcoming analysis of the A defor- 
mation under the classical perspective and also in the 
newer perspective [89j |. 

It is clear that the study of the A form factors is a very 
interesting topic, even if experimentally arduous. The 
only experimental source of information about the A de- 
formation comes today from the 7-ZV — > A quadrupole 
transition. However, this reaction depends on orbital an- 
gular momentum components in both the nucleon and of 
the A [6^, H3] • It is then very important to have an inde- 
pendent estimate of the angular momentum components 
in the A (leading to a deformation) [zl, IzE H2 ■ This is 
the main motive for our focus in this paper on the direct 
7A — > A reaction. 

In our previous work a first covariant spectator con- 
stituent quark model was developed for the nucleon, the 
7 A -> A transition and the A [j, |3i| [al HI]. In the 
covariant spectator quark model a baryon is described 
as a system of three constituent quarks, with a off-shell 
quark free to interact with a electromagnetic field, while 
the quark-pair (diquark) acts as a spectator on-shell par- 
ticle. Confinement is described effectively, with a quark- 
diquark vertex assumed to have a zero at the singularity 
of the three-quark propagator. The quarks are dressed, 
having an electromagnetic form factor with a behavior 
consistent with vector meson dominance. Within this 
first model we calculated the A form factors considering 
the A wave function parameterized simply by an S-wave 
state for the quark-diquark system. In this case only 
Geo and Gmi have nonzero contributions because of the 
symmetry assumed for the system. In the work reported 
here we extend the previous model to the inclusion of 
the D-states, following Ref. Q. It becomes then possible 
to predict, without further parameter fitting, nonzero re- 
sults also for Gei and Gm3- Although our results regard 
only the valence quark contributions, they compare well 
with a variety of lattice QCD data, besides the available 
experimental data. 

We organize this paper in the following main sections: 
In Sec. |TI]we introduce the basic definitions used to cal- 
culate the A electromagnetic form factors. In Sec. IIIII wc 
review the covariant spectator quark model and present 
the A wave function. In Sec. lIVI we derive the A form fac- 
tors within our model. In Sec. [V] we show the numerical 
results. Finally, in Sec. IVII we draw some conclusions. 



II. A ELECTROMAGNETIC FORM FACTORS 

The interaction of a photon with an on-mass-shell A 
isobar, with initial four-momentum P_ and final four- 
momentum P + , can be parameterized in terms of the 



3 



electromagnetic current |90j [9J 



J" 



-w a (P+) 



+ 



2^V 



4M A 



2M A 



7 p 

M p -)> 
(i) 



where Ma is the A mass, w a is the Rarita-Schwinger 
spin 3/2 state, F*{0) = e A , with e A = 0, ±l,+2 the 
four possible A charge values. The four functions F* 
(i = 1,...,4) of Q 2 define the four A form factors [in 
appendix^ we discuss the general form of (JTJ] ■ They can 
be combined to form the two electric and two magnetic 
multipoles, which have a direct physical interpretation: 
the electric charge (E0) and quadrupole (E2), and the 
magnetic dipole (Ml) and octupole (M3) form factors, 
defined as 0, [H |j 



Geo(Q 2 ) = [FZ-tF;] (1 + -t 



[F*-tF:]t{1 + t) 



Gmi(Q 2 ) = [F? + FZ](1 + ~T 



[f; + f:}t(i + t) 



(2) 



(3) 



G E2 (Q 2 ) = [F* — tF%] 



._[F*-rfJ](l + r) (4) 



Gm3(Q 2 ) = [F?+FZ] 



-[F 3 *+F 4 *](l + r) 



(5) 



where a and b are respectively the admixture coefficients 
for the two D-states: D3 (quark core spin 3/2) and Dl 
(quark core spin 1/2). N is the normalization constant. 
To characterize V^a we took the covariant spectator the- 
ory for a quark and a on-mass-shell diquark system, as 
proposed in Ref. [4(. 

From its symmetry properties, the J = 3/2 quark- 
diquark S-state (quark core S = 3/2 and orbital angular 
momentum L = 0) has the general form 

* s (P,fc) = ~MP,k) (T-e*) £p w a (P). (9) 

as proposed in Refs. 0, HI)]- I n this equation ips is a 
scalar function describing the momentum distribution of 
the quark-diquark system in terms of the A and the di- 
quark moment, P and k respectively; e* P is the polar- 
ization state associated with the diquark spin j92( (the 
polarization index was omitted in this short-hand nota- 
tion) and wp is the Rarita-Schwinger state. The isospin- 
| to isospin-i transition operators are represented by 
Ti (i = x,y,z) 0,1311101, and ^*(I Z ) is the spin-1 
diquark state (a combination of u and d quarks). The 
isospin operator (T • acts on the A isospin states. 
(See Ref. [HJ for details.) The normalization of the wave 
function (j9|) to unity implies for the scalar wave function 
the normalization condition 



Jk 



(10) 



where P = (Ma, 0,0,0) and the covariant integral over 
the momentum k is defined in Eq. (1251) below. 

The A D-states in the spectator formalism were intro- 
duced and explained in detail in Ref. [4j . Their definition 
is made through the 2?-state operator defined as 



where the dimensionless factor 

Q 2 



4M| 



(6) 



was introduced. The static magnetic dipole (ha), electric 
quadrupole (Qa) and magnetic octupole (0a) moments 
are defined in the Q 2 = limit, as 



HA = 



2M A 



Gmi(0) 



Qa - J jtGe2(0) 



Oa 



2Mi 



G M3 (0). 



(7) 



III. DELTA WAVE FUNCTIONS 

In our model the A wave function is described by an 
mixture of an S-state and two D-state components for 
the quark-diquark system 



a* 



D3, 



(8) 



where 



V^{P ± ,k ± ) = mi-^~gf, 



p, . u 



k% = k 



M 



po 



(11) 



(12) 



is the covariant generalization of the three-momentum in 
A rest frame, and 



papP 



(13) 



generalizes the unity operator. The variables P_ (P+) 
describe the initial (final) momentum, and in the same 
way k^ (k + ) refer to the initial (final) state. 

Using the operator V a generic spin 3/2 D-state is 



W a (P,fc;A) =V a P{P,k)w[){P,\). 



(14) 



It decomposes into a sum of two D-states, associated with 
the two possible spin states of the quark core, 5 = 1/2 
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or S = 3/2 (core-spin states ). To separate the two core- 
spin states we consider the two projection projectors: 



where 



77 M£ ' 



(15) 
(16) 

(17) 



The properties of the core-spin projectors are known in 
the literature [95ll96lj. 

The amplitudes for the states associated to the core- 
spin S — 1/2 (Dl-state) and 5 = 3/2 (D3-state) are 
defined as 

^di(P, k) = —3ipDi(T ■ e*)eTpW Dla (P, k) (18) 
*d 3 (P, k) = -3ip m (T ■ e*)er P W D3a (P, k), (19) 

where 

W Dla (P,k;X A ) = (V 1/2 ) a0 W (P,k;\ A ) (20) 
W D3a (P,k;X A ) = (V 3/2 ) a(j W?(P,k;X A ). (21) 

As before, (T ■ acts on the A isospin states. For 
simplicity, in our notation the diquark polarization (A) 
and A spin projections (Aa) are suppressed in the wave 
functions. Also ipni = ^Di{P,k), ipu^ = ipD 3 (P,k) are 
scalar wave functions. The factor —3 was introduced for 
convenience. The normalization of the wave functions is 



&D2s(P,k)} 2 } 



1, 



(22) 



and ensures that the wave function ^J A satisfies the 
charge condition, 

I *A(P,fc)ii(0) 7 °* A (P,fc) = e A . (23) 

A Jk 

The variable P represents P = P+ = P- in the A rest 
frame (Q 2 = 0) and ji(0) — \ + \ T z is the (quark) charge 
operator at Q 2 — 0. See Ref. [J] for more details about 
the D-states. 



IV. CURRENT AND FORM FACTORS 

Within the covariant spectator theory, the electromag- 
netic current in impulse approximation, can be written 
as in Refs. 0, M, Hi 

J " = 3 E / 9A(P+,k)tf*±(P-,k), (24) 

\ J k 



It also includes the symmetrization of the quark current. 
Following Ref. ;83| the photon-quark interaction is rep- 
resented by 



Ji = 



77/1- 
6 



-fl-Ts)-' 1 



g J2+ + ^h-Tz 



2M 



N 



(26) 



where Mn is the nucleon mass and the isospin (t 3 ) and 
spin (7^) operators acts on the A initial and final states. 
The coefficients fc± with i = 1,2 are the quark form 
factors and function of Q 2 . The explicit form of fi± will 
be introduced later in Sec. [V] 

For numerical applications we define the isospin depen- 
dent functions e A and k A as 



~e A {Q 2 ) = \[h+{Q 2 ) + h-{Q 2 )f 3 ] 



h A {Q 2 



[h+{Q 2 



f2-(Q 2 )T 3 



M A 



(27) 
(28) 



and T 3 is the isospin- 1 matrix defined as 



T 3 = 3j2 T hsTi = 



- 3 








" 





1 














-1 





. 








-3 



(29) 



We can also write T3 as X3 = 2Mt, where Mt is the 
diagonal matrix with the A isospin projections: Mt = 
diag(+3/2,+l/2,-l/2,-3/2) j3l|. In the limit Q 2 = 
we use 

e A =\{l+fz), k a = l( K+ + K _f 3 )^. (30) 

The values of n± were fixed in Ref. [83[ by the nucleon 
magnetic moments. 

It is convenient also to define 



9A 



e A - tk a , 



/a = Sa + ka- 



(31) 



Similarly as for e A and ka we suppress the tilde for 
Q 2 = 0. 

As both the initial and final states depend of the di- 
quark polarization with the same polarization index, the 
product of both polarization vectors factors out in the 
current (|2"4")l . Then, since the initial and final state have 
the same mass (M A ) [831181192]. 



13* 



(32) 



becomes 



which sums in all intermediate diquark polarization and 
integrates over the diquark three-momentum 



d 3 k 



{2nf2E k 



(25) 



A a/3 



papf 



+2 



{p+ + p-) a {p+ + p-Y 



AMI 



Q 2 



(33) 
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Using Eq. (J5J we can decompose the transition current 
into the contributions from the different components of 
the wave function: 

A 

x 

+3aN 2 |^ 9 Da j$*s + *sj^D3^ 
+3bN 2 &Dij?*s + ^2 ^sj^Di I 



+N 2 [0(a 2 ) + 0{b 2 ) + 0{ab)] 



(34) 



where the last three terms are corrections from the tran- 
sitions between the same or different D-states. Those 
terms can be neglected if the admixture coefficients are 
small. We can write J M in a compact form 



(35) 



with explicit forms for the currents Jg, J^ 3 and J D \ 1 
presented in Appendix [Bl 

From the general form (JT|) and for the wave functions 
defined by Eqs. © and fI5 |) -([l9" ]l . we can write 



Geo(Q 2 ) 


= N 2 g A l s 




(36) 


Gmi(Q 2 ) 


= N 2 f A 


4 

as + -aim - 
5 


D 


(37) 


G E 2(Q 2 ) 


= 3(aN 2 )~g A ^ 

T 




(38) 


Gms(Q 2 ) 


= IaN 2 


a h 2b 

T T 




(39) 



where Is is the overlap between the initial and final S- 
statcs, 



Is 



(40) 



and the other overlap integrals, between the initial S and 
each of the final D-statcs, arc 



Idi 
Id3 



b(k + ,q+)^ D1 (P + ,k)Tjj s (P-,k) (41) 
b(k+, q + )^ m (P+,k)Tjj s (P-,k), (42) 



with 



ul ~ \ 3 (~k+-q+) 2 lr 2 
b(k + ,q + ) = -—^ -k + 



(43) 



The variable k + is defined by Eq. (TT2"1) ; the variable g+ 
is defined the same way (just by replacing k by q). The 
derivation of the previous expressions is presented in Ap- 
pendix [5] 



In Eqs. (fBl) ]) - (|3"9")l , we chose to have N 2 = 1 in order to 
reproduce the charge of the A. (the corrections due to 
the D-state to D-state transition are of the order of a 2 , 
b 2 and ab, as indicated by Eq. ([34)) ). 

The first observation to be made on the results (|36|) - 
(131)1) is that, to first order in a and b, the D-states do not 
contribute to Geo ■ The second observation is that the D- 
states provide the only nonvanishing contributions to the 
electric quadrupole and magnetic octupole form factors, 
Gei and Gm3 (opening the possibility for a satisfactory 
selection of models for the D-states, otherwise specially 
difficult due to the overshadowing of their effects in Gmi 
by the dominant S-waves). 

The definition of b(k + ,q + ) implies that 6(fc+,q+) = 
-k 2 Y 20 (k) when Q 2 = (see Ref. 0] for details). Since 
the wave functions are independent of z = cos 9 in that 
limit, the angular integration of b{k+ 1 g+) vanishes in the 
limit Q 2 = 0, and therefore the D-state overlap integrals 
vanish also. We stress that this result is model inde- 
pendent, since it is a consequence of the orthogonality 
between the S-state (L = 0) and the D-states (L = 2). 

Consider now the form factors Ge2 and Gm3- They are 
proportional to 1/r, which is infinite for Q~ = 0. But, as 
Id3 or Idi are themselves proportional to r, both form 
factors go to a finite value as r and Q 2 — > 0, ensuring 
no divergence. The proportionality Ie3^Ib\ ~ t ~ Q 2 
simply comes from the specific dependence of the scalar 
wave functions on the quark momentum. The proof of 
that proportionality is done in Appendix [Cj 

This is how a nonzero contribution from the D-state 
overlap integrals survives for each form factor, in spite of 
the orthogonality between S and D-states. This same be- 
havior (overlap integral I ~ Q 2 as Q 2 — > 0) was already 
observed also in the Coulomb quadrupole form factor of 
the 7^ — > A reaction Q . 

To conclude, the form factors Ge2 and Gm3 are a con- 
sequence of the D-states alone, and they depend on both 
the admixture coefficients, a and 6, and the particular 
momentum dependence of the D-wave components of the 
wave function. 

A remark about the predictions for the A case can 
also be made: In our approach Ge2(0) = and Q A <> = 0, 
because the D-state to D-state transitions were neglected. 
Once those currents are included we may expect a 
nonzero, although rather small, contribution to Ge2{®)- 

The formulas obtained for the form factors are con- 
siderably simplified in the case Q 2 = 0. Since in that 
limit the integrals associated to transitions to the A D- 
states vanish due to angular momentum, via 6(fc+, q+) ~ 
^20 (k), we are left with, setting A^ 2 — > 1 



Geo(0) = e A 

Gjifi(O) = /a 

G E 2(0) = Zae A I' m 

Gm3(0) = /a [alro 



(44) 
(45) 
(46) 
(47) 



2bl' m }. 

(The factor Is = 1 from the normalization of tps was 
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suppressed.) In these equations 



We then use HH: 



X' m = Urn 

r— >0 



T' D1 = lim 



1m{r) _ dl D3 
dr 
dip i 
dr 



T 

1di{t) 



(0) 
(0), 



where Ids and I^i are defined in Eqs. (|4"Tj) - (l4"2"j) . 

The several moments, corresponding to each of the 
four multipole form factors, are calculated from Eqs. (J7J. 
Note that since the magnetic moment is defined by 
Gmi(0) and Idi and Ids vanish at the origin, it fol- 
lows that the magnetic moment is fixed by the S-state 
alone. 



V. RESULTS 

The algebraic results for the form factors given in the 
previous section were obtained within the covariant spec- 
tator formalism, independently of the specific model in- 
gredients. In this section we start to specify the scalar 
wave functions and the quark current, required for an 
actual numerical application. 

For the quark current we took the parameterization 
inspired on vector meson dominance (VMD) pi 183. |84| : 



h±(Q 2 ) 
h±{Q 2 ) 



= A + 



c±Q 2 /Mj 

(l-rf±) 

1 + Q 2 /m 2 + 1 + QZ/Ml 



(1-A) 
+ Q 2 /ml 

d± 



(48) 



where m v represents a light vector meson m v — m p (or 
m u ), and Mh = 2Mm is the mass of an heavy vector me- 
son which simulates shorter range effects. The parameter 
A was adjusted to give the correct quark density number 
in deep inelastic scattering [H,[84[, A = 1.22. The coef- 
ficients c and d were fixed in a previous work to describe 
the nucleon elastic form factor data. We have c+ = 4.16, 
c_ = 1.16 and d + = d- = -0.686 [&|. The particu- 
lar combination (|4"5|) with c + ^ c_ gives f\+ ^ /i_ and 
therefore breaks isospin symmetry for Q 2 ^ 0, essential 
for a good description of the nucleon data. The anoma- 
lous magnetic moments n + = 1.639 and re_ = 1.823 gen- 
erate the nucleon magnetic moment exactly (83 ] . The pa- 
rameterization (|48|) was applied to the nucleon [H, [§3[ , 
A electromagnetic form factors 0, [lij as well as to the 
three 7 TV ->• A transition form factors (1 I8ll85l. 193] . 

The scalar wave functions depend on the kinematic and 
dimensionless variable x> as introduced in Refs. (3.13.183. 

lallaS: 



X = 



(M A - m s f - (P - kf 



4>s(P,k) 
^ D1 (P,k) = 



N s 



m s (ai + x) 3 
N D3 
nr*(a 2 + x) 4 

N m r 1 



3 S l(a 3 + x) 4 K + x) 4 



(50) 
(51) 
.(52) 



MAm s 



(49) 



The functions ([5U]) - (T5S]) were chosen for a good descrip- 
tion of the 77V -> A data 0,[11[8|]. The coefficient X D1 
in ippi was determined by imposing the condition that 
the A Dl-state and the nucleon S-state are orthogonal 

In the applications we consider two different parame- 
terizations for the A wave function, both previously de- 
rived in the spectator formalism. By comparing the re- 
sults of two models we can illustrate the sensitivity of 
the results to model building. The two parameteriza- 
tions are equivalent with respect to the dominant S-state, 
but differ substantially in the D-state contributions. In 
both models the D-states give minor contributions for 
the quadrupole 7./V — ¥ A transition form factors, since 
they are basically dominated by the pion cloud contribu- 
tions at low Q 2 , which can be parameterized by a simple 
analytical structure [ij . The first model applied here, la- 
beled Spectator 1 (Sp 1), was derived in Ref. (3 (and 
named model 4 there) to describe the jN — > A transi- 
tions in the physical region. The D-states were calibrated 
by fitting the experimental data. The second parameter- 
ization, labeled Spectator 2 (Sp 2), was given in Ref. (85j . 
Although it has the same structure as Sp 1, it was ad- 
justed to the lattice QCD data [98] for 7iV — > A, in a 
pion mass region where the pion cloud effects are negli- 
gible, and afterwards extended to the physical pion mass 
point. The results for the physical region are therefore 
independent of the pion cloud mechanisms. This fea- 
ture gives the latter model some robustness, making the 
procedure, in our opinion, more reliable, and indirectly, 
constrained by QCD as well. In the first model (Sp 1) 
there is an admixture of 0.88% for the D3-state and a 
larger admixture of 4.36% for the Dl-state. In the Sp 2 
model one obtains an admixture of 0.72% for both D3 
and Dl. These results supply us with the first quantita- 
tive idea of the sensitivity of the D-states to the fitting 
procedure. The Dl-state, specially, shows to be strongly 
sensitive to that procedure. 



A. Electric charge form factor 

The results for G E a(Q 2 ) are shown in Fig. [T] The re- 
sults for Sp 1 and Sp 2 are exactly the same because both 
models have the same S-state parameterization, and we 
chose N = 1. In the figure, for the A + case we com- 
pare our results with the physical extrapolation of the 
A + form factor from quenched lattice QCD data [73l99l|. 
For the other charge cases, in the absence of lattice data 
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one can use exact isospin symmetry, which amounts to 
take 



G E + + (Q 2 ) = 2G£ + (Q 2 ) = -2G£~(Q 2 ), 



(53) 



where for G ^n we used the dipole approximation of the 
lattice data |73l . l99j . Although in quenched approxima- 
tion this symmetry holds, there is no rigorous theoretical 
justification for it, since breaking of the isospin symmetry 
can be expected for Q 2 ^ 0, and was seen to be consis- 
tent with the nucleon data. In fact, the spectator models 
accommodate isospin asymmetry (through /i + =/= /i_ 
and /2+ 7^ h-)- This is why, in Fig. [21 where 



is compared directly with G^q , the differences are more 
significant than what can be observed in Fig. Q] for A + . 
As expected, the disagreement between — G^ and G^ 
increases with Q 2 . Finally, in Fig. E] we compare the re- 
sults for A with the neutron electric form factor. Note 
that the slope near Q 2 = is very similar, for both mod- 
els and Geu- 

Comparing our calculations with the lattice data (ex- 
trapolated to the physical region) from Ref. [73|, |99| we 



note, as observed previously in Ref. Q where only S- 
states were taken, that both spectator models differ from 
the lattice data for low Q 2 , but are significantly close to 
the lattice data at larger Q 2 . Nevertheless, our models 
cannot be simulated by a pure dipole dependence as the 
lattice QCD data can, and one notices that they have a 
slightly slower falloff with Q 2 , implying a smaller charge 
radius or, equivalently, a stronger charge concentration, 
than suggested by the lattice simulations. In our calcu- 
lations the squared charge radius of the A + is 0.29 fm 2 
in both models, to be compared with 0.48 fm 2 from the 
lattice simulation 01 . 



Electric charge radius: comparison with other models 



The analysis of the charge distribution is naturally 
done first in terms of the charge radius. For the charged 
A states, the expansion of the charge form factor in pow- 
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Quenched Wilson [74, 751 






0.425±0.011 






Dynamical Wilson [74, 751 






0.373±0.021 






Hybrid [74, 7^ 






0.411±0.028 






Quenched Eg] 




0.410±0.057 


0.410±0.057 





0.410±0.057 



TABLE I: Summary of existing theoretical and lattice results for < r% > (fm 2 )- For the nucleon one has r 2 = 0.766 fm 2 and 
r 2 = —0.116 fm 2 [52]. The proton electric radius in the Spectator models is r 2 , = 0.79 fm 2 (see model II from Ref. [83]). 
a Extrapolation to the physical point. 
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ers of Q 2 , (e A = Gbo(O) + 0) 
Geo(Q 2 ) = G EO (0) 



i-^<4o>+o(Q 4 ) 



(54) 



defines the charge squared radius as 

2 6 dGso 



G BO (0) dQ 2 



(55) 
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For the case of neutral states (A ), Geo(0) = 0, and we 
use 



< r 2 > _ a dG EQ 
<r ^>- 6 dQ 2 



(56) 



Q 2 =0 



In some works Eq. (I56[) [with no normalization to Geq(0) 
included] instead of Eq. (|55l) . defines the charge radius 
also for charged particles. The definition (1551) has the 
advantage of being suitable to higher order form factors, 
namely Gmi, without loss of generality. For instance, for 
the proton, the electric charge and the magnetic dipole 
form factors have the same dipole dependence on Q 2 at 
low Q 2 , however with no normalization at Q 2 — as in 
Eq. (|56|) . one is led to very different electric charge and 
magnetic dipole radii. In addition, with Eq. (|55j) the A - 
and the A ++ radii can be directly compared with A + . In 
particular, in a model with exact isospin symmetry [see 
Eq. ([53] the charge radius is equal for all the charged A 
states if Eq. ([55]) is used. 

The results from our two models, together with a sum- 
mary of the literature, is presented in Table |TJ For a 
better interpretation of the results we write in the first 
column the prediction of each calculation for the squared 
charge radius of the proton, when it is available. Our 
prediction for the A + radius is 0.29 fin . This number is 
below the value predicted by a variety of models. 

The first estimate of the A charge radius was based on 
NRQM and suggested r\ + ~ 1.06 fm 2 [H[, implying a 
larger spatial distribution for the A than for the proton 
(r 2 ~ 0.77 fm 2 ). That effect was traditionally explained 
as a result of the repulsive hyperfine interaction of the 
quarks in spin-triplet state, in contrast with the attrac- 
tive effect in the singlet state 
at the time of this explanation, 
evidence of that fact. 



281 ] . However, there was, 
no direct experimental 



The A charge ra dius was also computed within the 
MIT bag model [HEo^. One should mention for com- 
pleteness other approaches, as a Skyrme model [HI], a 
quark-soliton Model (QSM) a field theory quark 

model (FT QM) [55(], a relativistic quark model (RQM) 
[56| . a constituent quark model (CQM) (68|, a gen- 
eral parameterization of QCD combined with large-iV c 
(GP/Large N c ) [58] a Goldstone boson exchange (GBE) 
model [ 431] , a chiral soliton model (xQSM) [13] and a chi- 
ral perturbation theory model (xPT) [49( . Although the 
results differ from model to model, they all have in com- 
mon the feature that the A has a larger charge radius 
than the pro ton. The exception is the xPT approach 

fromRef.gl- 

The same trend is seen in CQMs. By taking a nonrel- 
ativistic CQM with two-body exchanges, Buchmann et 
al. [31J obtains for the A + electric radius r\ + squared, 



t 2 
r A+ 



(57) 



where r 2 and r 2 are the proton and neutron squared ra- 
dius. In the same model r^ = 0. Using recent values 
for the nucleon radii (r 2 = 0.77 fm 2 ; r 2 = -0.116 fm 2 ) 
0.88 fm 2 , which, in any case, is 



[52j | one obtains r 2 A+ 
larger than r 2 . 



The relation (|57|) was improved by Dillon and Mor- 
purgo using a general parameterization (GP) of QCD (5?| 
with corrections of the order of 10%-20%. Ameliorating 
upon Ref. [3l|, Ref- HH establishes that the contribu- 
tions of the impulse approximation are dominant to the 
charge radius, and that the two-body currents associated 
with the quark-antiquark pairs are only a small correc- 
tion (« 0.03 fm 2 ). Using an operator method based on 
a GP [57J and the large- N c limit, Buchamnn and Lebed 
[58| related the four A charge radii with the neutron and 



10 



proton radii. Those results and their uncertainty are also 



shown in Table [U In this case ; 



A similar model, 



' A+ ~ ' p m 

as the GBE model [43[, predicts however a smaller A 
charge radius {r\ + ~ 0.7 fm 2 ). Finally, a calculation 



within xQSM suggests r 



A + 



0.8 fm 2 m 



Interestingly, the most recent estimates of the A elec- 



tric squared radius [31|, |33j, |42 



|62[ suggest r A ps 0.8 



fm , slightly larger, but not significantly larger than the 
experimental result for the proton r 2 ~ 0.77 fm 2 . In 
general we can say that the several calculations in the 
literature predict r A+ > r 2 . 

Our prediction for the A + radius (0.29 fm 2 ) is below 
all the model calculations identified above. In general 
these last ones also overpredict the lattice calculations 
values (~ 0.40 fm 2 and ~ 0.48 fm 2 from 0,1^). Only 
the estimate of xPT [2^] is close to these lattice results. 

As for the A charge radius, it is interesting to note 
that it is very close to the neutron radius (r 2 = —0.116 
fm 2 ) in our two models. This comes as no surprise, 
since, as noticed in Fig. G^o and the neutron elec- 
tric form factor, have almost the same shape for low Q 2 . 
The result is also consistent with partial quenched xPT 

(PQ X PT) m. 



2. Electric charge radius: comparison with the lattice QCD 

results 



Interestingly, our results are comparable in magnitude 
with the lattice QCD results in Refs. [73476|. The lattice 
results shown in Table U correspond to the lowest pion 
mass m n taken by each method (m^ ~ 300, 400 MeV). 

It is expected that the lattice results for the charge ra- 
dius increase as m n approaches the physical pion mass 
value, since it is what happens for the proton, where r 2 
diverges as m w — > 0. In fact, the A + lattice data [7a - l76l | 
follows approximately the behavior of the prot on charge 
radius as parameterized by xPT in Ref. jl0l| . This is 
seen in the left panel of Fig. 2J The vertical line shown 
indicates the inelastic threshold for pion production, de- 
fined by m-n- = Ma — Mm- When this inelastic channel is 
crossed from above we can conjecture the suppression of 
A + charge radius below that point, as suggested by the 
A magnetic moment studies within chiral effective field 
theory 0, [zl, Hl[ . In the graph we include also the result 



of xQSM in Ref. [62j at the physical point, r 2 A+ = 0.794 
fm 2 , which is close to the proton charge radius. 

In the right panel of the Fig. |4] we illustrate that the 
A + lattice charge radius results [76|, although compara- 
ble with the proton charge radius lattice data |l02| . are 
slightly larger, for the same lattice QCD conditions (lat- 
tice space a = 0.128 fm) and approximation (quenched). 

Some care must be taken in comparing model results 
like ours, and lattice QCD results. First, there are sev- 
eral methods to simulate the quark field in lattice, corre- 
sponding to different choices for the a ction s , as the Wil- 
son (quenched or dynamical) [zjlzl [Tol |lo||, Clover 
[Zl, Il02j or hybrid [zl |7]| actions. The different choices 



lead to different results. Besides, lattice methods can 
differ in the calibration needed to make the connection 
with the physical limit. That can be done by adjus ting 
the lattice spacing using t he S ommer method [l02j . or 
the physical nucleon mass 1041 ] . The two methods lead 
to very d ifferent predictions for the proton charge radius 

To finish the discussion on the charge radius, we men- 
tion that the method from Ref. (4||, based on quenched 
and partial quenched chiral perturbation theory, is useful 
to extrapolate lattice QCD calculations to the physical 
limit, and predicts results similar to the ones from our 
models Sp 1 and Sp 2. 



3. Overview of the electric charge form factor 

Apart from some deviation at low Q 2 , our models give 
a reasonable description of the lattice data in all their 
range, when lattice uncertainties are considered. 

Our two models (Sp 1 and Sp 2) gives the same result, 
because there are no D-state contributions for Geo and 
we have used TV 2 = 1. 

In comparison to other quark models, ours underes- 
timate the charge radius, while it is very close to the 
lattice data from Refs. [73M751 l99l |. We predict, as the 
lattice data does, a stronger charge concentration at the 
origin for the A than for the proton. Still, we obtain a 
slightly weaker charge concentration than lattice QCD. 
The quality of the agreement with the quenched lattice 
data increases for high Q 2 , where the meson cloud con- 
tribution should be smaller. This supports the idea that 
at least the valence sector is well described by our co- 
variant spectator model. However, some corrections are 
expected from the pion cloud at low Q 2 , and they are 
not yet included in our valence quark model. In princi- 
ple, they will increase the charge radius. Another effect 
to be added in the future is the inclusion of the D-state 
to D-state transition. But this effect is expected to be 
small since it is proportional to the D-state to D-state 
contributions of the order of a 2 , b 2 and ab, with a and b 
already very small. 



B. Magnetic dipole form factor 

The magnetic dipole form factor results are shown in 
Fig. [5j for the four charged A states. Only the magnetic 
moments of A ++ and A + are experimentally known, 
and their data points are represented in the grap hs for 
Q 2 = 0, according to ^a = Ga/^O)^— [Si— llOl] . The 
extrapolation of the quenched lattice Qdb data to the 
physical point exists only for the A + case. In that case 
the dependence on Q 2 of our results is directly compared 
to that extrapolation. For the other charge cases, in or- 
der to compare with the lattice QCD data we use the 
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Dynamical Wilson [74, 751 




2.35±0.16 






Hybrid [74, 75J 




3.05±0.24 






Quenched [76] 


5.28±0.92 


2.64±0.46 





-2.64±0.46 


Background [80] b 


4.85±0.16 


3.15±0.08 


0.0013±0.0210 


-2.42±0.08 



TABLE II: Summary of recent results for Gmi(O) compared with the experimental data and the SU(6) result. A compilation 
of earlier results can be found in Ref. To obtain )ia in /ijv units (nucleon magneton) use /ia = Gmi (0) ^^fiN- 
a Extrapolation for the physical point: b Ref. [8^] uses fi A - — ~^^++- 

The lattice results from Refs. [7414761 , |80| ] include no extrapolation to the physical point. The lattice results correspond 
respectively to m n = 411 MeV for quenched Wilson, — 384 MeV for the dynamical Wilson and m w — 353 MeV for the 
Hybrid action JzHzl]; rrw = 306 MeV for quenched from Ref. JzE] and = 366 MeV from Ref. H3|. 



relation based on exact isospin symmetry 

Gtn(Q 2 ) = 2GtA(g 2 ) - -2Gtn(Q 2 ), (58) 

together with the dipole approximation for the (Q 2 ) 
lattice data. The result extracted from the lattice data 
[78j using xPT [8l| for Q 2 — is also shown in the graphs. 

Our predictions for A + are very close to the lattice 
data. In detail, our results are 7% above the lattice data 
for low Q 2 . However, the lattice errorbars are larger than 
the ones for the Geo calculation, and our results are close 
to the limit of the errorbars. For high Q 2 the agreement 
is very good. In the limit Q 2 = our results are con- 
sistent with both the experimental points and with the 
lattice data extrapolation from Refs. [7^, HH, [99[. As 
found before in the calculation with S-waves only Q, 
the exception to this good agreement occurs in the A 



case, where [8l| predicts a sign different from ours. This 
disagreement is not problematic, since in that work the 
theoretical uncertainty was not explicitly evaluated, and 
Ha. is expected to be smaller relatively to the charged 
cases. 

In Fig. [6] we compare the effect of the D-states to Gmi, 
for both models Sp 1 and Sp 2. We conclude that the 
effects of the D-states are small although they increase 
with Q 2 . The small magnitude of these effects is easily 
understood from Eq. (|37[) with small admixture coeffi- 
cients. 



1. Magnetic moment: comparison with other models 

Our predictions for the magnetic moments are /xa++ = 
5.08/iat and fi&+ — 2.49/ia? for both Sp 1 and Sp 2. The 
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FIG. 5: Gmi form factor. Quenched lattice data from Ref. [73l . |99H . For A+ we include the experimental result from Ref. [f|: 
Gmi(0) = 3.54^|'|g(stat+sys) ± 3.94(theor) Q. The theoretical uncertainty for the A+ case is not presented in the figure. 



reason why both models lead to exactly the same result 
is that Gmi{0) is determined only by the S-state param- 
eterization, which is the same for the two models. There 
are some discrepancies between the different experiments, 
and the Particle Data Group @ reports the interval of 
3.7-7.5 for the A ++ magnetic moment. Our results are in 
good agreement with that result. Also the result for A + 
is consistent if the isospin symmetry is used to estimate 
the experimental value of fi&+ . 

The magnetic moment of A + is traditionally compared 
with the proton magnetic moment fi p . In the heavy 
quark limit (static) approximation they coincide. The 
SU(6) prediction Gjvfi(O) = 3.65 [IB] implies that, when 
expressed in the nuclcon magneton units (/xjv = 2 m n )' 



the A + and the proton have the same magnetic moment 
[ /iA+ = 3.65 _e_ = 3. 65 _e_ = 2.79^}. Our results 

are Gmi(0) = 3.29 in both models. 



for the A H 

In a previous work [3j, we compared our results 
for Gmi(0) without D- waves with several formalisms, 
namely relativistic quark models |29| : QCD sum rules 
[32I [34 |. chiral and soliton models [30, El| and dynam- 



ical reaction models (DM) based on hadron degrees of 
freedom [O, Now we compare our calculations, with 
D-waves included, with a comprehensive update of the 
more recent results based on xQSM [13], large- N c [13], 
SU(2) xPTjp, X PT [H, large iV c - X PT PQ X PT 
[45| , GBE [43| , QCD sum rules [69| , a quark models with 
sea quark contributions , a hypercentral quark model 
(HCQM) El, a Gauge/String Duality model [Hj], a re- 
cent DM 42 j and the U-spin symmetry 46]. For com- 
pleteness, to this list we add the CQM with one and two- 
body currents [3l|. The predictions for Gmi{0) from all 
these works, together with our own results, are displayed 
in Table HT1 The experimental and the SU(6) results are 
also included. Finally, we present the lattice QCD results 

[zi4zi, Lii, LzE HU ■ 

In general, all the different models are consistent with 
the available experimental information. 

The PDG interval for Gmi(O) corresponds to A ++ and 
is 4.8 - 9.8. One has 2.4 - 4.9 for A + assuming the 
isospin symmetry. The symmetric interval hold for A - , 
in the same approximation. Only Refs. 0, E3] for A~ 
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Quenched Wilson [74, _75] 






0.240±0.023 






Dynamical Wilson [74, .75] 






0.230±0.017 






Hybrid [74, 75] 






0.250±0.033 






Quenched [76] 




0.410±0.057 


0.410±0.057 





0.410±0.057 



TABLE III: Summary of existing theoretical and lattice results for < r\ n > (fm 2 ). Ref. [7(| the electric and magnetic radius 
are assumed to be the same. For the nucleon we can consider the result associated with the dipole form A 2 = 0.71 GeV 2 . The 
proton magnetic radius in the Spectator models is r\ lv = 0.74 fm 2 as extracted from Ref. [83l ] (model II). 
a Extrapolation for the physical point. 



are clearly out of the interval. 

As for the A case, SU(6) and exact isospin symme- 
try predicts Gmi(0) = 0, and our models, and in general 
all models predict contributions an order of magnitude 
smaller than the magnetic moment of the ch arg ed A or 
even zero (excluding the results from DM [42|], SU(2) 
XPT 47] xPT [49] and large iV c -xPT When ^ A o 

is not zero a negative value is frequently obtained, with 
a few exceptions. Note that the measurement of ^a 
is technically extremely difficult, since on top of the very 
short A lifetime, there is in addition the problem of track- 
ing neutral particles. For all these reasons our knowledge 
of had will probably be restricted in the near future to 
values from theoretical models and lattice simulations. 



2. Magnetic moment: comparison with the lattice QCD 
results 



Since the lattice simulations of the A magnetic mo- 
ment are performed for heavy quark masses, to compare 
with phenomenological models the extrapolation to the 
physical point is needed. 

Fortunately, there are nowadays methods based on 
xPT that can be used for such extrapolations [3, [zl, [U ■ 
The magnetic moment from R ef. I8ll i s the result of a chi- 
ral extrapolation. References [73l . lT9j j use a simple func- 
tional dependence on the pion mass to extrapolate to the 
physical point. All the other lattice results in Table |TT] 
refer to pions masses larger than the physical one. 



The dependence of the A + magnetic moment on the 
pion mass is presented in Fig. [7] in addition to the re- 
sult from models Sp 1 and Sp 2 at the physical point. 
The lattice data in that figure comprises the quenched 
results from [74T - [7ol l78l . |79| , and the unquenched results 
from [74l |75| obtained with the Wilson and hybrid ac- 
tions [8 0| | using Clover fermions and also the background 
field method [73, [ZSl, H3|- The chiral extrapolations of 



Cloet [81| using the pion loop corrections, and Lee [79( 
using a simple analytical form, is also on the graph. The 
inelastic point (m^ < Ma — Mn) is represented by the 
vertical line. The solid line in the figure corresponds only 
to the analytical part of the chiral extrapolation from 
Cloet et al. [8l|, but the analytical line and the full re- 
sult lead to similar results for fi&+ a l the physical point 
[8~fl | with an accuracy better that 10%. The extrapolation 
given by Cloet et al. compares well with our results. 

The main conclusion from the figure is that there are 
large discrepancies between several lattice calculations. 
This makes harder to draw conclusions from the compar- 
ison of our results with the available lattice data. The 
unquenched results from Aubin et al. [80] for the A + ex- 
ceed all the other calculations, including the analytical 
contribution to the chiral extrapolation [81]. Similarly 
to Ref. [79j, in the Aubin results a significant deviation 
of the isospin symmetry (see Eq. (JSSJ) is manifest in the 
A ++ and A + predictions. Note however that the Aubin 
results for while different from the A + results, 

are similar to the results from other calculations. The sig- 
nificant violation of the condition h^A++ ~ Ma+j shows 
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FIG. 6: Comparing the S-state approach with the effects of the 
D-states in model Sp 1 and Sp 2. Lattice data from Ref. ff^ . 



that the Aubin calculation violates strongly isospin sym- 
metry, given by Eq. (l58l) . 

To understand the discrepancies between the differ- 
ent lattice calculations, it is important to realize that 
besides the extrapolation to the physical point, certain 
lattice calculation of Gmi demand an additional extrap- 
olation of Q 2 down to the Q 2 = point, given that the 
minimum lattice momentum Q 2 is non-zero. The most 
common method uses agiven analytical form for the ex- 
trapolation. A dipole 73(, or an exponential falloff form 
[74L [75l | have been used. Anot her possibility, used by 
the Adelaide group [zH [ll, Il02j | , applies the scaling be- 
tween the electric charge form factor and the magnetic 

GmiIQ 2 ) _ G E0 (Q 2 



dipole form factor, 



This assumption 



Gjvii(O) Geo(0) 

implies that the magnetic dipole radius squared is the 
same as the electric charge radius. Note, however that, 
although the scaling approximation is valid for the pro- 
ton form factors [83[ , there is no reason to believe that it 
holds also for the A system. 

As a complete different alternative, the background 
field method [z3, [7^, H3] takes the three-quark system 
interacting at rest (Q 2 — 0) with a static magnetic field. 



Considering the energy shift in the system, which is pro- 
portional to the magnetic dipole and also with the mag- 
netic field, the magnetic moment is determined avoiding 
the uncertainty present in the extrapolation of the form 
factor method. The different results of [§(J may be a con- 
sequence of the combination of a strong magnetic field 
with a not so large lattice volume used in the numeric 
lattice simulations. 

The other problem to consider when comparing be- 
tween and to different lattice calculations lies in estimat- 
ing how the quenched approximation deviates from the 
results from the full QCD calculation, with the same pion 
mass. 

In the SU(6) limit (exact SU(3) flavor symmetry) gives 
MA+ = /%>• This limit is obtained in quenched QCD 
when and light and strange quarks have the same mass, 
i.e. rriK = m^, with uik the kaon mass. In a recent 
quenched simulation [76| the SU(3) symmetry point cor- 
responds to ~ 700 MeV (m 2 = 0.485 GeV 2 ) [76j]. 
Below this point (m^ < 700 MeV) [i v > ha+ ■ This result 
can be understood analyzing the pion cloud contributions 
in quenched and in full QCD. As explained in Ref. [76| . 
considering again A + , for < 700 MeV, those con- 
tributions are negative in quenched approximation while 
positive in full QCD. 

Therefore, for A + in the region < 700 MeV, the 
quenched result underestimates not only full QCD, but 
even the the core contribution, which is already below 
the full QCD result. This is an artifact o f the quenched 
approximation also observed in Refs. [711 Il05l|. 

As for /ia°! the extrapolations of Cloet [8l| and Lee 
[79l ] give different signs. The Aubin result for = 366 
MeV [8(| is positive but the errorbars are so large that 
it can also be consistent with zero and negative values. 
Note that the results of Lee [79| are quenched and the 
results of Aubin [80j are unquenched. 



3. Magnetic dipole radius 



We define the magnetic dipole radius squared < 



A/1 



> 



by means of the expansion of Gmi specified analogously 
to Eq. (f5"5"f. The results of the magnetic dipole radius 
squared are listed on Table IIIII for different models. At 
low Q 2 they are similar to the ones obtained for the elec- 
tric charge radius, suggesting a scaling between the elec- 
tric charge and magnetic dipole form factors, as it was 
already observed for the proton form factors. 

In the comparison between our results and others, the 
main difference is seen in the A magnetic dipole radius, 
where we obtain a value larger than the radii predicted 
by all models discussed till this point. This can be just a 
consequence of the small value of Gmi for Q 2 = 0. As it 
happened for the electric radius, the results for the mag- 
netic dipole radius of our two models, as well as the lat- 
tice QCD data, underestimate the results obtained with 
other quark models (see results of Refs. [H, |H, |62| ) ■ This 
can be a consequence of not including the pion cloud, ex- 
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FIG. 7: A + magnetic moment from lattice calculations. 
Quenched lattice data (labeled Q) from [fUzl [?|, [ll H3 . Un- 
quenched lattice data from 74, 75, SOj] - The solid line is the ana- 
lytical contribution for fJ.^+ as derived in Ref. [8lll . The result of 
Sp 1 and Sp 2 (filled >) is also represented. 



pected to be important at low Q 2 



4- Overview of the magnetic dipole form factor 

Our results indicate that the S-state contributions are 
dominant for the A magnetic dipole form factors. In the 
previous section we saw that the same holds for the elec- 
tric charge form factor. The D-states contribute only 
with small corrections. In spite of pion cloud effects 
not being included explicitly here, we predict magnetic 
dipole distributions similar to the lattice results. When 
compared with other quark models we predict a larger 
concentration of the magnetic dipole distribution at the 
origin. 



C. Direct comparison with lattice QCD 

We have compared our results for the magnetic dipole 
form factor with the extrapolation to the physical pion 
mass point of the quenched lattice QCD data. But 
quenched data include pion cloud effects only partially, 
and lattice results based on quenched approximation are 
known to underestimate the full QCD result, as dis- 
cussed. Fortunately, all A form factors were recently 
also evaluated with unquenched methods. Those include 
sea quark effects from up and down quarks (Wilson ac- 
tion), and also from the strange quark (hybrid action) 
[74l L75[ . For those calculations there are no extrapola- 
tions to the physical pion mass limit available. In these 
cases we could only compare our models directly to the 
published lattice data. Figure [5] shows the results for 
Geo and Gmi corresponding to the lowest pion mass in 
each case: 411 MeV (quenched), 384 MeV (Wilson ac- 
tion) and 353 MeV (hybrid action). In addition to the 
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FIG. 8: Lattice data from Refs. [73 
other results. 



75, 76 



underestimates the 



results of Refs. [75j we include the point corresponding 
to = 306 MeV from Ref. [76j|. 

Surprisingly, our two models also reproduce the mag- 
nitude of the lattice QCD data even in the nonphysical 
pion mass region. The agreement between our predic- 
tions for Gmi and the lattice data at low Q 2 explains 



why the < 



Ml 



> values obtained with our models are 



close to the lattice values, as seen in Table Unl Note that 
although different methods give very similar results for 
Gmi they exceed slightly the physical point (one data 
point from Ref. [76j is the exception, probably as a con- 
sequence of the small pion mass considered). 



D. Electric quadrupole form factor 

Because of Eqs. (|3^)) - (j57| . the charge and magnetic 
dipole form factors are essentially determined by the S- 
state component of the A wave function. But without 
D-states in the A wave function the electric quadrupole 
form factor is identically zero. As such, Ge2 measures 
a transition between the S-state and the D3-state (S — 
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TABLE IV: Summary of existing theoretical and lattice results for Ge2(0). To obtain the electric quadrupole moment use 
Qa = G E 2(0)-rfj- with M A = 6.24 fin" 1 . 
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Quenched Wilson [75J 




0.336+0.096 






Dynamical Wilson [75J 




0.134+0.106 






Hybrid [75] 











TABLE V: Summary of results for < r| 2 > (fm 2 ). 



3/2, L = 2) of the A [see Eq. ([55])]. Differently said, 
the D3-state induces a deformation in the A system that 
builds up a nonzero electric quadrupole proportional to 
the D3-state admixture parameter a. 



1. Electric quadrupole moment 

The electrical quadrupole moment Qa is denned as 
Qa = Ge2{0)j§t- The shape of the A can be inter- 
preted according to the sign of Qa, if one relates in the 
Breit frame the charge distribution to the form factors in 
the nonrelativistic limit. For a positive charge, a defor- 
mation extended along the equatorial region (Qa < 0) 
corresponds to an oblate distribution (pancakelike) and 
a deformation along the polar axis (Qa > 0) to a pro- 
late distribution (cigarlike). Alternative interpretations 
of deformation can be introduced, and give different in- 



sight on the structure [zl, [7|| . More details about this 
issue of deformation are presented in a separate work [89j j . 

Our results for Ge2 at Q 2 = are shown in Table IIVI 
Since Gei depends only on the state D3, Ge2(0) is de- 
termined by I' m which is —7.00 (Sp 1) or —6.65 (Sp 2). 
Besides, the D3-state admixture coefficents are not very 
different (0.88% versus 0.72%) and therefore both mod- 
els give similar values for Ge2(0)- In our case, Qa° = 
because the D-state to D-state contributions were ne- 
glected. But nonzero although small contribution for 
Qao may appear otherwise. 

In Table HVl we make comparisons to other works, start- 
ing with nonrelativistic quark models (NRQM) e.g. the 
classic Isgur model [63]. In those models the quark D- 
states are a consequence of the hyperfine interaction in 
the presence of a confinement mechanism. Then Qa 
was estimated in terms of the admixture coefficients for 
the symmetric, antisymmetric and with mixed symmetry 
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FIG. 9: Form factor G E2 - The A+ lattice QCD points correspond to m n = 353 MeV from [7j|. The NRQM (©) is the result of Eq. J59j 
and [|3l; x?T (o) is from The xPT result for A is out of the scale (see Table ITV)| . 



components in the S and D-state wave functions, together 
with a confinement parameter associated to the spherical 
harmonic oscillator (interpreted as the quark core radius) 
[III [3l| , and assuming, as we do here, that the transition 
between D-states is negligible. 

In this type of NRQM description, where the valence 
quark degrees of freedom and the electromagnetic inter- 
action are reduced to the impulse approximation (which 
includes only the one-body current), it was possible to 
relate the A cha rge d istribution with the neutron charge 
distribution [66], Il06| according to 



Jimp) 



9 2 



(59) 



where r 2 is the neutron squared radius, and the index 
(imp) holds for impulse approximation (one-body cur- 
rent). Equation (|59p is a parameter free relation, since 
it is independent of the admixture coefficients and of 
the confinement parameter. Considering a recent esti- 
mate r 2 = -0.116 fm 2 [13, one has Q A+ ~ -0.0464 

fm 2 , or (0) — —1-81, in close agreement to the re- 
sult of Ref. [66| . Similar numerical results were obtained 



by Buchmann et al. [3l| with a different combination of 
values for the admixture coefficients together with the 
confining parameter. 

To estimate the nonvalence contribution for Qa absent 
in the formulas above, Buchmann et al. [3l| included 
two-body currents for the quark-antiquark production 
mechanisms, and obtained 



(60) 



In this case no D-state admixture is considered. Al- 
though it is a consequence of the constituent quark for- 
malism, this result has again the feature of being pa- 
rameter independent |3llj. and could even be obtained in 
the large N c limit [601 ] - Later, Eq. (1601) was improved 
by means of a GP of QCD [53, [53, [61J, for the inclu- 
sion of higher order terms. From the 7^V — > A electric 
quadrupole data [H| the value G B2 (0) = -7.02 ± 4.05 
was obtained. Note that this result, and already the re- 
sult from Eq. (IBTfl) [3~0 ], are large when compared to the 
other calculations compiled in Table HVl 

The other results in the literature are based on xPT 
[H [H, [ll , on a quark-meson coupling model (QMCM) 
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FIG. 10: Comparing Ge2 with lattice data from Refs [74. l75j| . Left panel: quenched data (open o); Wilson action (filled □) and hybrid 
action (filled o). The quenched result from [76H for m-jr = 306 MeV is also included (filled o). Right panel: interpolation of the lattice data 
[ZllZl and xQSM [fjl]; NRQM (o) [H and QCDSR (v) [H results are also shown. 



65|. on a Skyrme model [67J. and more recently xQSM 
62| and QCD sum rules [63] ■ For completeness we 
mention as well calculations of PQxPT in the heavy 
quark mass limit |45l . l7l| , and relations between the elec- 
tric quadrupole moments for different charge A states 

Mm 

From lattice QCD there are earlier calculations from 
Leinweber [78| , and recent simulations [62|, l74r{76j based 
on different approaches: quenched, dynamical Wilson (u 
and d sea quarks) and hybrid action (also s quarks in- 
cluded). The lattice results are also presented in Table 

El 

Our model compares well with the NRQM from 
Refs. (HI, HH HU based only on the valence quark de- 
grees of freedom (see Table IIV[) , although it overshoots 
the QCD sum rules [69| . Mixed descriptions, involving 
quarks and mesons, in particular pions, as degrees of free- 
dom as xPT jHSEa. QMCM [g| and xQSM H suggest 
that the meson cloud can be important, although not as 
significant as Eq. ([t)0"f seems to imply. 

2. Dependence of Ge2 with Q 2 

Before we compare our results with lattice QCD data 
we investigate the dependence of Ge2 on Q 2 . The results 
for Geo{Q 2 ) are presented in Fig. |9l for the models Sp l 
and Sp 2, for all the A charge cases. For A + the results 
are also compared with the more complete unquenched 
lattice QCD simulation (hybrid action) for the lightest 
pion mass available (m, = 353 MeV) [75| . As a reference 
result we include also the xPT calculation from [49| and 
the NRQM from Eq. [The xPT prediction for A 

is out of the scale although the errorbars are consistent 
with our result (zero)]. The hybrid action considers u 
and d quarks as valence quarks, and u, d and s as sea 
quarks. Therefore pion cloud contributions are included, 



although the results are not extrapolated to the physical 
pion mass. 

In the left panel of Fig. [TU]our models for A + electric 
quadrupole are compared with other lattice QCD meth- 
ods, namely the quenched calculation with m T = 41 1 
MeV and the dynamical Wilson action calculation with 
TOtt = 384 MeV, and again the hybrid action results 
for m n = 353 MeV. All the cases shown correspond to 
the smallest pion mass values of each method. In the 
same graph, also the lattice analytical extrapolations for 
Q 2 = assuming an exponential dependence in Q 2 are 
included [Note the significant errorbars]. 

In the right panel of Fig. [TU] we compare Sp 1 and 
Sp 2 with the dipole approximation of the model xQSM 
from [62], that provides results for Ge2 as function of 
Q 2 . Additionally, the figure depicts the results from the 
analytical approximation of lattice QCD of Refs. [zHlzIi; 
as well as the NRQM and QCDSR [69| results for Q 2 = 0. 

From Fig. [TU]we conclude that our predictions are close 
to the hybrid action and the dynamical Wilson action 
results. Still, our results are larger, in absolute value, 
than the quenched data, particularly in the low Q 2 region 
(note that the dynamical hybrid data are consistent with 
the Wilson and quenched data, within their statistical 
errors). The difference between our results and different 
lattice approaches seem to suggest that contrarily to the 
leading form factors Geo and Gmi 0], the meson cloud 
effect may be important for Ge2 ■ Also, the model %QSM, 
where the contribution of the meson cloud (sea quark 
contributions) near Q 2 = is around 70% 62] falls faster 
that those models as Q 2 increases. With the exception 
of Ref. 73] where an extrapolation to the physical point 
was published, there are no extrapolations of Ge2 results 
for the physical point. 

In the left panel of Fig. [10] we depict also the lattice 
data point of Ref. f7(| estimated in quenched approxi- 
mation for m v = 306 MeV at Q 2 = 0.230 GeV 2 . In the 
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limit of the errorbars, this point is also consistent with 
both our models (Sp 1 and Sp 2), 

To finish the discussion about Ge2(Q 2 ) we present in 
TablefVlthe result of our models and the ones from lattice 
QCD for < 

r E2 >• Al so f° r this observable Sp 1 and Sp 
2 are similar, as we would expect from the graphs from 
Fig. IH1 In both models the radius is in agreement in 
magnitude with the results from lattice. 



3. Overview of the electric quadrupole form factor 

The A electric quadrupole form factor is a measure of 
the deformation in the charge distribution, and is pro- 
portional to the transition coefficient between S and D3- 
states. Except for the lattice QCD calculations [73l - l75j |. 
and for the model xQSM [62|, previous studies of the de- 
pendence of the A electric quadrupole form factors on Q 2 
are almost nonexistent. There are however predictions of 
several models at the photon point. Our model is con- 
sistent with most of these calculations based on valence 
quark degrees of freedom. 

Our results agree also with the dynamical lattice QCD 
simulations with (Wilson and hybrid actions) for the 
lightest pion cases. Lattice QCD data seems to suggest 
that sea quark effects play a more important role for Gei 
than for Geo and Gmi, but the statistics must still be 
increased to allow more definitive conclusions. 

The similarity of the D3-state parameterization in both 
our models Sp 1 and Sp 2 leads to very close results for 
the electric quadrupole form factor. We conclude that 
the quadrupole electric form factor is not suitable to dis- 
criminate between our models. 



E. Magnetic octupole form factor 

The third order form factor, Gm3, did not attract as 
much attention as the first order (Geo, Gmi) and the 
second order (Ge2) form factors. Even for the magnetic 
octupole moment, Oa, there is no experimental informa- 
tion whatsoever, although Oa was estimated considering 
the earlier NRQM [28[. Interestingly, though, the lattice 
QCD simulations in Ref. [73| raised the question whether 
it is possible to determine the sign of Oa accurately, or 
equivalently, with respect to the magnetic dipole distri- 
bution, whether the A is prolate (Oa positive) or oblate 
(Oa negative). The first lattice QCD simulation 78 1 
was not accurate enough to answer this question. The 
recent work of the MIT-Nicosia group [73[ based on the 
quenched approximation increased the hope that Ga/3 
or at least its sign can determined unequivocally by lat- 
tice QCD. This was a challenge for both MIT-Nicosia 
[75[ and Adelaide groups [76j]. The results from Ref. (7(| 
(which are quenched) suggest a negative sign (i.e. a pro- 
late distribution) for Gm3- Nevertheless, one has to check 
whether the results for Gm3 are not significantly affected 
by the sign of the pion loop contributions, as it happens 



one was Z^, 3 = 
I' D3 — —6.65, I'rji 



for Gmi (see discussion in Sec. IV B[) . One should note 
that those results correspond to large mass pions, and no 
physical extrapolation was performed yet. 

From the theoretical point of view, the publication of 
lattice QCD results for Gm3 triggered the application of 
quark models to compute Oa and Gm3- Buchmann and 
Henley [2J using a CQM with pion cloud, two-body cur- 
rents and the GP formalism [57| predicted Oa+ = —0.012 
efm 3 . The magnetic octupole was also recently predicted 
by QCD sum rules [69[ and xPT [Zi|. Predictions for 
Gm3 (Q 2 ), based on the xSQM are expected for a near 
future [62;]. xPT predicts a nonzero octupole moment 
only at two meson loops level [3(J. Within PQxPT in 
the heavy quark mass limit at one meson loop level, one 
has Ga/3(0) = [45ll7~Lll. consistently with the quenched 
QCD data from Ref. [76[ for large pion mass values. 

A summary of the predictions for Gm3(0), including 
our results from models Sp 1 and Sp 2, is presented in 
Table ED 

In the covariant spectator quark model, according to 
Eq. (j4"?| , G 4/3(0) depends only on the integrals T' m , I' Dl 
and the admixture coefficients a and b. For model Sp 1 
7.00, T m = 1.59. As for model Sp 2, 
0.24. For the model Sp 1 the large 
contribution of the Dl-state (4%) cancels the D3 contri- 
bution leading to an almost zero value to Gm3(0). This 
entails a dramatic difference between the results of mod- 
els Sp 1 and Sp 2, as we can see from Table PVTl Model Sp 
1 gives negligible contributions for Gm3(0). The result 
from Sp 2 lies between the negligible predictions of QCD 
sum rules, and the large value of Buchmann l2f based on 
a pion cloud model and the GP formalism 0, [53, [6(| • For 
completeness, and following the previous sections, the re- 
sults for < r M3 > are also presented in Table IVlTl 

The differences between models Sp 1 and Sp 2 for Gm3 
are clear from Fig. [TTJ For Q 2 ^ the cancellation be- 
tween Dl and D3 contributions in model Sp 1 still exists, 
although it is not as spectacular as in the Q 2 = case. 
We note also the change of sign in Sp 2 for Q 2 w 0.4 
GeV 2 . It would be interesting in the future to compare 
the Q 2 dependence of model Sp 2 with alternative mod- 
els or with lattice data. The estimates of QCD sum rules 
Hi, XPT [H and GP 2] for Q 2 = are also included 
in the graphs. For the A + we show the lattice data of 
Ref. [76j for the pion mass m n = 306 MeV at Q 2 = 0.230 
GeV 2 . In the graph we represent also the range of varia- 
tion obtained in Ref. [73[ from calculations with several 
pion masses. 

The results for A + from model Sp 2 underestimate in 
absolute value the lattice QCD data point from Ref. [76[ 
for m-K = 306 MeV, but it is consistent with the central 
value of xPT [49| . We note, although, that the error- 
bands are significant. The predictions of Sp 2 are also 
consistent the data from Ref. (73|, within an interval as- 
sociated with several pion masses. It will be interesting 
in the future to verify if the quenched approximation af- 
fects crucially Gm3 and if the extrapolation to the phys- 
ical point introduces significant corrections. To achieve 
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FIG. 11: G M 3 form factor. The results for QCDSR (V) and xPT (filledo) [U for Q 2 = are shown. Lattice data for A+ from 
Refs. [ll,[75| (m n = 306 MeV) for Q 2 = 0.230 GeV 2 , and Alexandrou et al. [73| for Q 2 = 0.42 GeV 2 . The error bar for the data from 
Alexandrou et al. represents the interval of values associated with several pion masses. For A+ the result of GP is also included. 



G A f3(0) 


A++ 


A+ 


A 


A~ 


QCDSR [69] 


-0.0925 ± 0.0234 


-0.0462 ± 0.0117 





0.0462 ± 0.0117 


GP [2] 


-11.68 


-5.84 





5.84 


X PT [49] 


-4.07±9.49 


-2.03±4.76 


0.059± 0.036 


2.26±4.52 


Spectator 1 


-0.049 


-0.024 


0.00084 


0.026 


Spectator 2 


-3.51 


-1.72 


0.064 


1.85 



TABLE VI: Summary of existing theoretical results for Gm3(0). To obtain the magnetic octupole moment use Oa = 
GM^Wjfpr, with M A = 6.24 far 1 . 



that goal the errorbands must be reduced. At that point, 
Sp 2 can then be better tested. 

We can conclude that, contrarily to the form factors 
analyzed before (electric charge, magnetic dipole and 
electric quadrupole) the magnetic octupole form factor is 
extremely sensitive to the D-state parameterization. As 
for our calculation enabling a decision in favor of Sp 1 or 
Sp 2, we consider model Sp 2 our best model, because 



it is better constrained, i.e. it is constrained by not only 
the experimental but also the lattice data of the reaction 
7^ A [85| (see discussion at the beginning of the sec- 
tion). Model Sp 2 predicts a much stronger deformation 
from the symmetric magnetic dipole distribution. It is 
also peculiar that the Sp 2 predictions for Gm3(0) are 
similar to the predictions of Ge2(0)- This implies an al- 
most identical deformation for the electric charge and for 
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< r 2 M3 > 


A++ 


A+ 


A 


A" 


Spectator 1 


~150 


~150 


~150 


~150 


Spectator 2 


1.2 


1.1 


2.0 


1.2 



TABLE VII: Summary of results < r 2 M3 > (fm 2 ). 



the magnetic dipole distribution, if one uses the classical 
perspective of deformation. 

A look at the magnetic octupole charge radius shows 
also a dramatic difference between the two models. 
Model Sp 2 predicts typical (although slightly large) val- 
ues, but model Sp 1 breaks all the scales. That result 
might however be due to the equation used for < rf vI3 >, 
which is the analogue of Eq. (|55[) . and breaks down in 
the limit when Ga/3(0) gets close to zero. But it might 
also be a consequence of the pion cloud effects not being 
explicit in this work. 



1. Overview of the magnetic octupole form factor 

Model Sp 2, which is our best model, gives a result for 
the magnetic octupole closer to the available lattice QCD 
data. We admit that model Sp 2 is not the last word 
in the A wave function parameterization. Future lat- 
tice QCD simulations with increasingly smaller system- 
atic errors may still discard the parameterization of Sp 
2 85] . Note that the quadrupole transition form factors 
can nowadays be experimentally measured for that reac- 
tion (the interaction time in 7^ — > A is larger than in 
the reaction 7 A — » A). A direct test of the parameteriza- 
tion will be an even more accurate study of the 7-/V — > A 
reaction, once more precise data is provided. By using 
quadrupole form factors (electric and Coulomb) at the 
physical point [ 4 | . 18411 or, in alternative, their data from 
lattice QCD [74|, Ha], the valence contributions of the A 
wave function, and in particular the small Dl and D3- 
state contributions, can eventually be better constrained, 
to further refine our the predictions for the Ge2 and Gm3 
A form factors. 



VI. CONCLUSIONS 

In this work we extended the spectator quark model 
to include the A wave function D-states into the matrix 
elements of the electromagnetic current, in first order in 
the admixture coefficients of those states. This allowed 
us to evaluate all the four A elastic form factors. We took 
two models pre viously calibrated by the two nucleon form 
factor data 83j and by the form factors for the 7^ — > A 
reaction [3, 85]. The two models differ in the way they 
describe the jN — > A quadrupole data. The first model 
(Sp 1) was derived by fitting the experimental data to 
the sum of the valence quark contribution with the pion 



cloud contribution [4]. The second model (Sp 2) was ob- 
tained by adjusting only the valence quark contributions 
to that reaction to the lattice QCD |85j |. data while still 
describing the physical data. Because in the last case 
we fixed the A D-state parameters independently of the 
pion cloud mechanisms, in our view model Sp 2 is more 
reliable and robust than Sp 1. 

Our predictions, based on the valence quark degrees 
of freedom alone, reproduce well the main features of 
the lattice results for Geo and Gmi when extrapolated 
to the physical point, except for a some deviation at low 
Q 2 . In particular, the agreement for high Q 2 (say Q 2 > 1 
GeV 2 ) is excellent, indicating that pion cloud effects (not 
included explicitly in our model) are small for high Q 2 
values. This confirms that the valence quarks are the 
dominant effect in the Geo and Gmi form factors, as 
anticipated in a previous work Q without the explicit 
contribution of the D-states. In the low Q 2 region, how- 
ever, our results for the A suggest a smaller charge and 
magnetic charge radii, when compared with the proton. 
This can be a limitation caused by the absence of the chi- 
ral behavior (or pion cloud effects). Alternatively, it can 
also represent the real physical case if, as for the A mag- 
netic moment, there is a suppression of the charge radius 
due to the nonanalytical contribution of the inelastic re- 
gion (m^ < Ma — Mj\r). Further lattice QCD studies, 
preferably unquenched studies, are needed to clarify this 
point further. 

The inclusion of the D-states in the A wave function 
is required to explain the electrical quadrupole and mag- 
netic octupole form factors. The D3-state (spin 3/2 state 
with an orbital D-wave) state is needed for a nonvanish- 
ing electric quadrupole moment for the A. The Dl-state 
(spin 1/2 with an orbital D-wave), as well as the D3-state 
(spin 3/2 with an orbital D-wave), generates nonvanish- 
ing contributions for the magnetic octupole form factor. 
By comparing our results with the available variety of 
lattice results, we conclude that valence quark contribu- 
tions are of the same magnitude as the total result. This 
shows that to establish the magnitude of the contribu- 
tions of the pion cloud to these observables requires more 
progress, with more accurate quenched and un quen ched 
lattice data. In its absence for the near future [l07j . we 
may estimate the pion cloud effects similarly to w hat was 
done for the nucleon and the octet in Ref. |l08j. That 



work concluded that, although significant for the nucleon 
magnetic moments, pion cloud effects give corrections of 
less than 8% for the quark anomalous magnetic moments. 

Another interesting result to be further explored in the 
future is that contrary to Ge2, Gm3 is very sensitive to 
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model parameterization. This oversensitivity can be ob- 
served in the Gm3 radius. The value of Gm3(0) depends 
strongly on the D-states admixture and on the overlap 
between S and Dl-states. We obtained an almost com- 
plete cancellation of Gm3 in a particular model (Sp 1). 
If Gm3 (0) is in fact very small, then the study of D-state 
to D-state transitions may became important and should 
also be explicitly evaluated. Since the direct measure- 
ment of a third order multipole is technically difficult, 
a way out is, in our opinion, the use of the jN — > A 
(experimental and lattice) data for an even more accu- 
rate determination of the D-state admixture following 
Refs. HH, to predict O a as done here. 
Our best model (Sp 2) gives for A + 



G M i(0)=3.27, 
G M3 (0) = -1.72. 



G B2 (0) = -1.71 



(61) 



These values correspond to the multipoles: 

Ha+ = 2.49/ijv, Q A + = -0.044 efm 2 

£> A+ = -0.0035 efm 3 . (62) 

As the lattice QCD data for Gm3 is still affected by a 
large uncertainty, our results for Oa+ are independent 
predictions to be tested in the future. 

The results (joT]) - (|o^|l . give an oblate deformation for 
the electric charge and magnetic dipole distributions, 
considering the classical interpretation of deformation. 
In Ref. [89( we compare this analysis with the one based 
on transverse density distributions. 

The dependence of the Gm3 form factor on Q 2 that we 
obtain is very atypical, and it will be interesting if future 
lattice QCD calculations or alternative hadronic models 
may confirm our results. It can also happen that forth- 
coming lattice calculations may require a reparametriza- 
tion of our model, which would affect in turn our results 
for the 7 A'' — > A (where the quadrupole form factors, 
electric and Coulomb, are small and sensitive to the D- 
states.) 

In the future, we can extend the model t o the lattice 
QCD regime as suggested in Refs. [H, HI, Il09j |. That 
possibility was already explored in previous applications 
to the nucleon [HI and 7iV ->■ A transition [II,!!!]. The 
covariant spectator quark model defines a tool that can 
be used to extrapolate the lattice results both to the 
Q 2 — > and to the physical pion mass limits. Another 
strength of this formalism is the possibility of its exten- 
sion to the strange quark sector, enabling the application 
of the model to the octet and decuplet baryon syst ems . 
In particular, the application to the il~ form factors [l09j | 
is very promising, since meson cloud effects (kaon cloud) 
are expected to be very small in that system. Our pre- 
dictions could then be compared with lattice simulations 
nowada ys p erformed at the strange quark physical mass 

[7ll^liTq. 
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Appendix A: Generic form to the current 

Consider the A electromagnetic current J M 

J" = WcV^Wf}. (Al) 

The most general form, consistent with parity, G-parity 
[94j . gauge invariance (current conservation) and the 
properties 

(P± -M A )w a {P ± ) = 

w a (P+)P% = 0, F%(P_)=0 

l a w a = 0, (A2) 



reads 



= (ag af) + bq a q f} ) ^ 
+ (cg af) + dq a qP)P» 



-e g 



(A3) 



The coefficients a, b, c, d, e and / are general functions 
of Q 2 . Time reversal invariance implies that r(P + ,P_) 
must be symmetric in the permutation P + ^ P_ (or 
q — > — q), and then we conclude that 

/ = -e. (A4) 

Using Eq. (|A4|) . one has 

= (ag a(S + bq a qP) + 
(cg a(S + dq a q f} )P» + 

+ e ( 5 <V - q a g^) ■ (A5) 

To obtain the final form we use the identity due to Fear- 
ing (gj: 



g^(P + +P. 



^Mi + Q 2 ^ 



2 Ma 



-g~"T 



a a a /3 
w -4- _ q y ; 



M A 



and the Gordon decomposition 



(P++P-Y _ „ io* v q v 
2 Ma 7 2M A ' 



(A6) 



(A7) 
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which holds between Rarita-Schwinger states and solu- 
tions of the Dirac equation. Finally, one has 



F*g^ + F, 



F*g ali + F 4 



q a q fi 
3 4Mj 
q a q f3 
SMI 



r 



l<J^q u 



2 Ma 

where F* are linear combinations of a. b, c and e. 



(A8) 



Because the S and D3-states have the same core-spin, 
in general 

A Jk 
A Jk 

although the result becomes zero for Q 2 = 0. 



Appendix B: A form factors 



1. Current Jg 



The current Jg was written in Ref. [3J using the A 
S-state wave function of 1841: 



a. Transition S — > D3 



where 



G S E0 (Q 2 ) = g A Zs 
Gmi{Q 2 ) = /a^s 
G| 2 (<9 2 ) = G S M3 {Q Z 



Jk 



(Bl) 
(B2) 
(B3) 



(B4) 



As Gg2 = Gm3 = there is no electric quadrupole nei- 
ther magnetic octupole originated by the A S-states. 



2. Current J£ 3 

Consider Eqs. (|9|) and (fT9)l in a compact notation 

*s = -e?w a 1> s (T't U ) (B5) 
*D3 = -^*W ma ^ m {T ■ e*)- (B6) 

The transition between a A S and a A D3 state can be 
decomposed in two processes: 

• S in the initial state (current Jg D3 ) 

• D3 in the initial state (current J^ 3 s ). 

In the following we will use the properties of the spin 
S-states (|A2[) (S-state) and also 

P±W mp (P ± ) = M A W mp (P ± ) 
PlW D3 p(P±) = 0, T P W D 3i}(P±) = 0. (B7) 

Note that the last identity holds for spin states based 
in a spin 3/2 core dll|. Simple consequences of those 
relations are 



W D 3a 4wp = W a ^W D 3/3 = 0. 



(B8) 



The transition between an initial A S-state (momen- 
tum P_) and a final A D3-state (P+) can be written as 

J S,D3 = 3E / *B3(P+,k)j?* S (P-,k), (B9) 

where the spin and polarization indices are suppressed 
for simplicity. We can write (|B9[) as 



Js,D3 = 3j{[W ma A»w p ] ife^g) A^, (BIO) 
where ip£ 3 = tp m (P + ,k), ipg = ip s (P-,k) and 



eA7 + ka- 



2M A 



(Bll) 



In the last equation we use the definitions of §a and ka 
from Eqs. (|2T ]) -([2"8" j) . 

Using the Gordon decomposition (|A7jl . we write the 
current (|B10I) as 



(P++P-) 



-K s ,m, (B12) 



2Af A 

where /a is defined by Eq. (|3"Tj) . and 

Js,D3 = 3 J { [W D3 a -fw p ] iP+^s } A^ (B13) 
n s .D3 = i [ {[W D3a w p ]iP+ 3 iPs}A a P. (B14) 

Jk 

The expression for A a ° is given by Eq. (|3"3"|) . 

To reduce the current Jg D3 to a form close to the stan- 
dard one (HJ), we work the spin algebra using the prop- 
erties (|A2|) and (|B7[) . We also use the covariant result of 
the integration in the azimuthal variable tp [4j: 



Vfy+^g = R^{q,P + ) / b(k+,q+)i>£ 3 il>s> 

k Jk 

(B15) 



where 



b(k + ,q^ 



3 (fc+ ■ g +) 2 ljr 2 
2 # 



2 fc ; ( B16 ) 



R {q , P+ )- — -- g+ . 



(B17) 
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The variable <?+ follows Eq. and q+ comes from 

Eq. (fT2|) with k replaced by q. In the limit Q 2 — (initial 
and final state in the A rest frame) b(k, q) is reduced to 
~k 2 Y 2Q (k). 

At the end, one has 

1 + 2t q a qP 



Ji 



S,D3 



1 + T Q 2 



.rf + g"?^ 



1 g<V P» 



1 + tM a * 



1 + t Q 2 M A 



S,D3 



3 + 2r gV , Q/3 
The common integral is defined by 



wp x2+ 3 . 



-^D3 = / KKA+)^m^s- 



(B18) 
Wfs x 2+ 3 , 
(B19) 

(B20) 



where 



One can prove now that: 



X D3 = I b (k-,Q-)^ DS f s - 



7"+ 7 — — 7" 



(B28) 



(B29) 



One starts with the two integrals with explicit arguments 
written in the Breit frame. Then we conclude that the 
integrand of 2^ 3 is a function of z and the integrand of 
is a function of — z. Replacing z — > —z we obtain 
identity in the Breit frame. As the integrals are covariant, 
the result holds for any frame. 



c. Final result 



Transition D3 



The transition between an initial A D3-state (momen- 
tum P_) and a final A S-state (P+) can be written as 

J D3,S = 3 E / *s(P + ,k)j?* D3 (P-,k). (B21) 

Once again, we suppress the spin and polarization indices 
for simplicity. We can write (|B21I) as 



J D3,S 



3 / {[w a A^W D 3p]^m} (B22) 



where 4>D3 = i'D3(P-,k), tpg = ips(P+:k). The opera- 
tor A 1 is given by Eq. (|B11[) . 

Using again the Gordon decomposition (|A7j) . we obtain 



J D3,S 



D3,S 



2M A 



-7tD3,S, 



(B23) 



where /a is defined by Eq. (PTj) . and 

■7m,s = s [ [M L W mp ] A^ (B24) 

Jk 

K m>s = 3 ( {{WacWm^+^g} A a ?. (B25) 



To reduce the current (|B23|) to a form close to the 
standard one (JT|), we work the spin algebra as before. 
We use again the covariant result of the integration in 
the azimuthal variable ip from Eq. (|B15|) with the obvious 
replacement of P+ — > P- (which implies k + — > fc_ , q + — >• 
q_ and g°f — > gl 13 )- One obtains 



^D3,S — 



1 + 2t q a qP 



1 + T 



1 V ^ + ffi^ 



1 



1 + t 9 M A 1 + t Q 2 M A 



TlD3,S = w Q 



3 + 2r q a q 



1 



Q 2 



(B26) 
wp x l m , 
wp x l m , (B27) 



Adding the two currents 



J _D3 ~~ J S,D3 J D3,S' 



we define two components (a and &) 
T M — 7 M + 7 M 

J D3 ~ J a i J b 1 



where 



D3 



7'' 
■ J b 



k a 71d3 



{P+ + P-Y 

2 Ma 



(B30) 

(B31) 

(B32) 
(B33) 



with 



D-3 



2, 



D3 



1 + 2T9V 
1+T Q 2 



1 1 

1 + t Ma 



7?.D3 — 22j 



, 1 q a q f3 {P + + P-Y 
3 + 2t gV 



D3 



1+T Q 



(B34) 
(B35) 



Considering the Fearing relation (IA6|) and the Gordon 
decomposition (|A7I) we can finally write 



TV 



D3 



1 + T 



ID3 X 



d(3 



3 q a q> 3 
' tAMI 
2 q a q !i 



t4M 2 a 



2M A 



(B36) 



From the previous relation and the definition of the form 
factors F*, we conclude that the form factors associated 
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with J£, from Eq. (|B32|) . arc 



Finally we take the sum of the a and b components. 









Using the definition for /a, we are left with 




F* a (Q 2 ) 


2 

, . fA^D3 
1 + T 


(B37) 


G D E l(Q 2 ) = 


(B53) 


F; a (Q 2 ) 


2 

-, . /a£d3 
1 + T 


(B38) 


5 


(B54) 




6 ? Id3 
= -, . /a 
1 + t r 


(B39) 


Gf 2 3 (Q 2 ) = Tgjf- 


(B55) 


F; a (Q 2 ) 


4 f Xl53 
= i . JA 

1 + T T 


(B40) 


GZ(Q 2 ) = h^f. 


(B56) 



To evaluate the multipole form factors we consider the 
transformations ©-©. First we note that F 2 * a = -F* a 

!7*a 3 J7i*e 

4 - ~2^3 



and FT a = -|K* Q . Then 



F* a - tK 



2 



2 T 



It follows 



^EO 


= 


(B41) 








s~ia 


= -t/a^D3 

5 


(B42) 


r~ia 


= 3/a 

r 


(B43) 


"M3 


? 2Tr>3 

= /A ■ 

r 


(B44) 



Take now the current J£ from Eq. (|B33[) . Using the 
Gordon decomposition (|A7|) and the form of TZd3 from 
Eq. (|B35|) one concludes that 



F? b (Q 2 ) = 2k a 1 D3 
f;\q 2 ) = -F* b {Q 2 ) 



(B45) 
(B46) 



F£ (Q ) = 2-—~ka — 

1 + T T 



(B47) 

7 

Ft(Q 2 ) = ~F?{Q 2 ). (B48) 

From the previous equations it is easy to conclude that 
p*b + p *b = Q) F *b + p *b = Q 



F* - TF* b = 2(1 + r)K A 2 m 

3 



Ff - T F? = 2(3 + 2t)~k a 1 D3 . 



And then, 



G b EO =0 
G MX = 

G^ 2 = -3(1 + t)k a — 
r 



G 



b _ 
M3 — 



0. 



(B49) 
(B50) 

(B51) 

(B52) 



Note that the electric form factors include the electric 
factor g A , while the magnetic form factors include the 
factor /a- To achieve this final form it was essential to 
include the terms associated with TZd3- the result of the 
overlap between the spin functions associated with core 
spins 3/2 (S and D3). Note that, although I03 vanishes 
for Q 2 = 0, £22. _> const as Q 2 0. 



3. Current 

Consider Eqs. © and (fT5j) in a compact notation 

*s = -e?Tu> ^sr(T • (B57) 
*£,i = -Ser^m^m^-e 1 *). (B58) 

The transition between a A S and a A Dl state can 
be decomposed into two processes: 

• S in the initial state (current J§ D1 ) 

• Dl in the initial state (current J c \ 1 s ). 

In the next calculations we use the properties of the 
spin S-states (|A2[) (S-state) and Dl-state Wd\0'- 

P±W D1I3 (P ± ) = M A W m p(P±) (B59) 

i£Wdip(P±)=0. (B60) 

See Ref. Q for details. 

a. Transition S — > Dl 

The transition between an initial A S-state (momen- 
tum P_) and a final A Dl-state (P+) can be written as 

Js,di=zY1 I *m{P+,k)3$* s {P-,k). (B61) 

A Jk 

We suppress the spin and polarization indicies for sim- 
plicity. We can write (|B61|) as 

j s,di = 3 / {[W ma A»wp] A ^ (B62) 



The magnetic contributions b vanishes for Q 2 = 0. 



where ^zn — ipDi(P+,k), ips = ipDi(P-,k) and is 
given by Eq. (|B11[) . 
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The Gordon decomposition (|A7|) of the current (|B62|) 
gives: 



J^m = 3/a / {[Wm^w^^s] ^ 

Jk 

~ xp+ + p -Y 



a/3 



(B63) 



" 3/ ^ 2 Ma 
where /a is defined by Eq. (f3"Tj) . and 

K = [ {[W Dla wp]iP m iPs}A a P. (B64) 

The orthogonality between the states Wdi a and wp gives 
ft = 0, (B65) 

equivalently to 



(B66) 



By working the spin algebra and using the covariant 
integration in the azimuthal variable <p [4|: 



Ik Jk 

(B67) 

with b and R defined respectively by (IB16I) and (|B17I) . 

At the end, we obtain 



J S.Dl 



1 



1+T 



fAl m 



^ + 2^4-7" - 2 



Q 2 



Q 2 Af A 



where 



X di = J KK^+^ix^s- 



(B68) 
wp . 

(B69) 



b. Transition Dl — > S 
We consider now the transition from Dl to S: 

Jdi,s = ^Y1 [ *s(P+,k)jf*Di(P-,k). (B70) 
We conclude that 

4,s = 3/aE / {[w a rWmp]i> m i>t}^ aP - (B71) 
a ^ 

where the orthogonality was used again. 
Working through the spin algebra, 



Ik Jk 

(B72) 



we are left with 



Wq ^ 5 Af a + Q 2 7 Q 2 A/a 



(B73) 



where 

J- x = / (B74) 
Similarly to Appendix IB 2 b| we can prove that 



7"+ 7 — = 7" 

X D1 — -^Dl — L D1- 



c. Final result 



(B75) 



Adding the two currents (jB68[) and (jB73|l 

J D1 — ,J, S,D1 "r J D1,S> 



we obtain 



J D1 



= h 



1 



(B76) 



(B77) 



+4 



7 M - 4- 

2 ' r>2 



2M A 



Wp 



Using the Fearing- Nozawa relation (|A6|) and the Gordon 
decomposition (|A7|) . we are left with 



1 



(B78) 



a/3 



afj + 2 



4Af A 



2 g<V 



r 4M A 



2A/ A 



wp 



Comparing the previous equation with Eq. (|24[) , we con- 
clude that 



F* = 


1 . IaIdit 

1+T 


(B79) 


f; = 


. IaIdi 

1 + T 


(B80) 


F* = 


2 /aZdi 

1+T 


(B81) 


p; = 


2 /a 

1+T T 


(B82) 



Note that, although F£ includes a factor 1/t = 
4A/ A /Q 2 , it is not singular. By definition of b(q,k) the 
integral Idi vanishes as Q 2 —> canceling the divergence 
in 1/t. The condition Idi = when Q 2 = is the result 
of the orthogonality between the states L = and L = 2. 
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Using the simple relations 

f*(q 2 )^tf;(q 2 ) 
f;(q 2 ) = tf:(q 2 ), 



(B83) 
(B84) 



we can write 



Finally 



F?(Q 2 )+FZ(Q 2 ) = -2f A l D1 

f*(q 2 )-tf;(q 2 ) = o 

F*(Q 2 ) + F 4 *(Q 2 ) = -4/ A ^i 

T 

F*(Q 2 )-rFUQ 2 )=0. 



2 ~ 

Gmi(Q ) = — c/a^di 
5 

G^{Q 3 ) = 

G M3(Q 2 ) = 

r 



(B85) 
(B86) 
(B87) 
(B88) 



In the limit Q 2 = only Gm3 is different from zero. 



4. All contributions 

Each of the four A form factors G a (a = 
E0,M1,E2,MS) is then the result of the three contri- 
butions calculated above: 



G a (Q 2 ) = N 2 [Gf (Q 2 ) + aG% 3 (Q 2 ) + bG™(Q 2 )] 



(B89) 



The final result becomes 



N J g A X 8 

r2, 



Geo{Q 2 

Gmi(Q 2 ) = N 2 f A l s + ^(aN 2 )f A 2 D3 -^(bN 2 )f A l D1 
HaN 2 )g, Xm 



G E 2(Q 2 ) 

G M3 (Q 2 ) = IaN 2 



1A- 

T 
%D3 



2b 



%D1 



(B90) 



Appendix C: Proof that To ~ Q 2 

Consider the overlap integral 



1 D = b{k + ,q+)^D{P+,k)^ S {P^ 1 k), (CI) 
Jk 

where D represents an arbitrary (D3 or Dl) D-state. The 
functions %pE and ?ps are respectively the D-state and 
the S-state scalar wave functions written in terms of the 
dimensionless variable \, as defined in Eq. (|49|) . That is 
the only assumption made in this appendix. The function 
b(k + ,q + ) is defined by Eq. (|4"5|l. 



Although the integral Id is covariant the analysis can 
be considerably simplified in a particular frame. We con- 
sider then the rest frame of the final baryon (mass M): 



P+ = (M, 0,0,0) 
P_ = (£,0,0,-|q|) 
g=(w,0,0,|q|), 



(C2) 



where 



E = y/M 2 + |q| 2 = 



2Af 2 + Q 2 
2M 



V(4M 2 + Q 2 )Q 2 _ EQ 



2M 
P+-q _ Q 2 
M 2M ' 



M 



(C3) 



In this frame we can write k + = (0, k x ,k y ,k z ) and q + = 
(0, 0, 0, |q|). In the final baryon rest frame k + and q+ has 
only spacial components. 

It is now trivial to conclude (k + ■ q + = — |q|fc z ) that 

b(k+, q+) = - k 2 ) = -k 2 Y 20 (z), (C4) 

where we used k z = \k\z, with z = cos 9. In the final 
baryon rest frame there is no Q 2 dependence in &, and 
only the scalar wave functions depend on Q 2 . 

The momentum dependence of the scalar wave appears 
through the functions P± ■ k, which become in the final 
baryon rest frame: 



P + -k = ME, 



(C5) 



EQ 



P_-k = EE s + \q\k z =EE s + -^kz. (06) 

In the last equation we used k = |k|. The previous equa- 
tions show that the angular dependence appears only in 
the S-state wave function. To perform the angular inte- 
gration (variable z) for we need only to consider the 
factor 



I 



dzY 20 (z)ip s (P-,k), 



(C7) 



where it is implicit that the integration is in the interval 
[-1,1]- 

To proceed, we need to use the momentum dependence 
of ipS- By taking the generic form 



MP-,k) 



Ns_ 
D n ' 



(C8) 



where n is a integer (n > 2), N$ is some normalization 
constant, and 

P-'k EQ 
D = a + 2—-2 = P + 2^ (09) 



where /3 = (a - 2) + 2 



EE 3 
Mm,. 
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Factorizing D according with 



EQ k 



(CIO) 



we can write 
11 

ao + fll [wp)* + aa [Wp 



z 2 + 



> (Gil) 



where 



a, = 2-? 



= 2 3 



(-n)(-n-l)..(-n-i + l) 



(C12) 

with j = 0, 1, ... For n = 3, for example: ao = 1, ai = —6, 
a2 = 24, etc. 

The angular integration (|C7j) is now considerably sim- 
plified 



dz 



{ Y2o(z) 



ao + 



(EQk\ (EQ 



\M 2 P 



\M 2 P 



) H } 



«2 



(EQ 
\M 2 



dzz 2 Y 2 o(z) 



0{E 4 Q 4 ), (C13) 



where the integral of the first term is zero [J Y2q{z) = 0], 
the second term vanishes because is the integral of an odd 
function (in a symmetric interval), and the same argu- 
ment holds for the forth them (coefficient 0,3). The next 



nonzero contribution appears only with the fifth term 
(coefficient 04). As J dz {z 2 Y2o(z)\ = the previous 
result can be written as 



T _ ±^2_ ( E k 

~ 15/3" \M 2 p 



Q 2 + 0(E 4 Q 4 ). 



(C14) 



We still have to perform the integration in k. However, 
because ~ ^^-m s = const as k— > 00, all integrals are 
well defined providing that the original integral is well 
defined. This means that integration in k preserves the 
decomposition presented above. 

For small Q 2 we need only to consider the leading de- 
pendence in Q 2 of each term. Using Po = (a — 2) + 2^-, 
one has 



4 a 2 ( k 



15 PS \MP 



Q 2 + 0(Q 4 )~Q 2 (C15) 



Therefore 



Id ~ Q 2 



(C16) 



which finishes the proof. 

Note that this derivation is independent of n. The 
only constraint to be imposed is that the power in -05 
is sufficient to assure the convergence of the integral (in 
k). An alternative derivation is possible in a different 
frame, using also the properties of ipD- For instance, in 
the initial baryon rest frame all angular dependence is 
in i))d- Considering ipjj with the same form of ips from 
Eq. (|C8j) . we obtain the same limit for small Q 2 using 
the procedure described before, as it should be. 
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